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We consider scalar-Gauss-Bonnet and modified Gauss-Bonnet gravities and reconstruct these 
theories from the universe expansion history. In particular, we are able to construct versions of 
those theories (with and without ordinary matter), in which the matter dominated era makes a 
transition to the cosmic acceleration epoch. In several of the cases under consideration, matter 
dominance and the deceleration-acceleration transition occur in the presence of matter only. The 
late-time acceleration epoch is described asymptotically by de Sitter space but may also correspond 
to an exact ACDM cosmology, having in both cases an effective equation of state parameter w close 
to —1. The one-loop effective action of modified Gauss-Bonnet gravity on the de Sitter background 
is evaluated and it is used to derive stability criteria for the ensuing de Sitter universe. 

PACS numbers: ll.25.-w, 95.36.-fx, 98.80.-k 



I. INTRODUCTION 



Modified gravity is a promising theory that has become a very attractive gravitational alternative for dark energy (for 
a recent review, see [HID, 9- It i s J to start, a powerful scheme. Indeed, depending on the specific model considered, 
modified gravity is able to realize any of the proposed scenarios that have been delimited by the observational 
constrains leading to cosmic acceleration: effective phantom models, cosmological constant theories or quintessence. 
Also, it can easily account for the different epochs in the evolution of the unverse, that start to emerge clearly from 
the observational data (for a recent review of these data and their comparison with dark energy models, see Q). The 
qualitative understanding of gravitational dark energy (see Q for a review) is quite simple: some gravitational terms 
different from the usual General Relativity ones may dominate at the very early or very late universe epochs. Thus, 
General Relativity seems to be only approximately valid, both at very early as well as at very late times. To have 
the possibility to explain — in a unified way as modified gravity effects — fundamental cosmological phenomena such 
as early time inflation and late time acceleration, is very appealing. Moreover, modified gravity has the possibility to 
solve the coincidence problem, too, and can also clarify the role of dark matter in the formation and evolution of the 
universe. 
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Among the different classes of modified gravities, a very interesting one is the family of the string-inspired gravities. 
String-inspired scalar-Gauss-Bonnet gravity has been suggested in Ref. || as a possibility for gravitational dark 
energy. The idea that these theories may lie close to effective models which may come from fundamental string theories 
is in itself very appealing. Some time ago, Gauss-Bonnet (GB) gravity was applied to the possible resolution of the 
initial singularity problem [6[. Another version of string-inspired gravity, namely the modified GB or F(G) theory 
can also play the role of gravitational dark energy. The investigation of different regimes of cosmic acceleration in 
such string-inspired gravity models has been carried out in refs. [auL H, H EH El, 13, EH • 

A very strong theoretical restriction on modified gravity is caused by the natural emergence of the known classical 
universe expansion history. In other words, as in the case of the usual ACDM cosmology, it must faithfully reproduce, 
to start with, the sequence of the well-established cosmological phases: radiation/matter dominance, deceleration- 
acceleration transition, and the cosmic speed-up. Recently, a reconstruction scheme for the case of f(R) gravity has 
been developed [H, EH E3] where such a cosmological sequence is seen to occur quite naturally for some models 
(for a review on reconstruction from the universe expansion history, see (TtI EI| ) . Having in mind the fundamental 
importance of the correct description of the past and the current universes, in this work we develop a reconstruction 
program for string-inspired gravity from the universe expansion history. Using the method developed in refs. fl3l . , 
it will be here demonstrated that scalar-Gauss-Bonnet or F(G) gravity can be reconstructed for any given FRW 
cosmology. Moreover, the role of string-inspired gravitational terms may be quite important, even in the matter 
dominated epoch. Examples of exact and/or approximate ACDM cosmology in the above modified gravity are 
discussed too. The stability of such form of FRW cosmology, which naturally develops an asymptotically de Sitter 
future is investigated. As de Sitter universe naturally occurs at the early or late times in such models, special attention 
is paid to the de Sitter universe. Using the results of the calculation of the one-loop effective action of F(G) gravity 
in the de Sitter space, the semi-classical stability of de Sitter space is investigated. 

The paper is organized as follows. In the next section we develop the reconstruction scenario for string-inspired, 
scalar-Gauss-Bonnet gravity. Elaborating on the approach of Refs. fl3l[l9l|. it is shown that the cosmological sequence 
of matter dominance, deceleration-acceleration transition and cosmic acceleration may occur in the specific version of 
that theory with some given potentials. It is also shown that an exact ACDM cosmology can indeed be reconstructed 
from such theory. Section three is devoted to the reconstruction of F(G) gravity from the universe expansion history. 
For the same two classes of FRW universes (an approximate and an exact one), ACDM cosmology is derived in 
versions of the theory where this cosmology naturally occurs. By introducing perturbations around such cosmological 
solutions, the stability of those two FRW universes is investigated at the classical level. It is shown that the very 
early universe may be unstable, which opens the very interesting possibility of a natural exit from inflation, while 
the late stage, e.g. the asymptotically de Sitter universe, can be stable. Sect. 4 is devoted to the study of the same 
reconstruction scenario for F(G) gravity in the presence of several types of usual matter. Two specific versions of 
F(G) gravity with matter are constructed where the matter dominance and deceleration-acceleration transition occur 
only with matter. After that, the accelerating universe becomes asymptotically (or even exactly) the de Sitter one. 

In Sect. 5 we discuss the following problem: does it exist any scenario which may still improve the above situation, 
with the emergence of the matter dominated era for scalar-Gauss-Bonnet or F(G) gravity? It is shown that such 
scenario may be realized: by adding the compensating dark energy which is relevant in the matter dominance epoch 
mainly, the cosmological sequence of matter dominance, deceleration-acceleration transition, and cosmic acceleration, 
do occur in the versions of such theories which, otherwise, do not contain the matter dominance. Some discussion 
and outlook are presented in the last section. 

Appendix A is devoted to the calculation of the one-loop effective action of F(G) gravity on a de Sitter background. 
Some technical remarks about the appearance of the multiplicative anomaly in those calculations are done. In 
Appendix B, this one-loop effective action is applied to study the semi-classical stability of de Sitter space in F(G) 
gravity. This is an alternative method to investigate stability, which is advantageous as compared with the very 
involved one that uses cosmological perturbations. Numerical calculations provide examples of modified Gauss-Bonnet 
gravity which have a stable de Sitter vacuum solution. 



II. RECONSTRUCTION OF THE SCALAR-GAUSS-BONNET THEORY FROM THE UNIVERSE 

EXPANSION HISTORY 



In the present section, we will show how the string-inspired, scalar-Gauss-Bonnet gravity theory, proposed as a 
dark energy model in Ref. can indeed be reconstructed, for any requested cosmology. The starting action is 



R 
2^2 



d^d^-V^-^G 



(1) 
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Here G is the Gauss-Bonnet invariant and the scalar field <j> is canonical in {l}. We now assume that the FRW 
universe with scale factor a(t): ds 2 = —dt 2 + a(t) 2 Y^,i=i —a (dx l ) and the scalar field 4> only depend on t : (f> = 4>(t). 
These assumptions are based on the observational data indicating that the spatial part of the universe is flat and the 
universe could be regarded to be homogeneous on large scales. In the early universe, the small inhomogeneity could 
have played an important role. For example, the large scale structure of the present universe could be generated from 
the very small initial fluctuation and the gravitational instability for the perturbation. In this paper, however, we are 
interested in the behavior of the universe on the scales larger than the 'large' scale structure and we now assume the 
homogeneity of the universe. Hence, we now drop the spatial coordinate (x l ) dependence of a and <p. In the following, 
the Hubble rate H is defined by a/ a. 

Then the FRW equations look like [EEl]: 



= -- 2 H 2 + U 2 + V{4>) + 2AH 3 ^Z y j ' , (2) 



/v 



2 dt 



dt 

o = ±(2H + 3tf 2 ) + \tf - v(<f>) - sH^j^ - ieura*y _ 16H Mm (3) 

and the scalar field equation is 

O = + 3fty + V'(0)+£(0)G. (4) 
Now G = 24 (tiH 2 + , and combining © and ©, one gets 

= \H + 2 - 8^^« - leHB^m + ^AlWl = ^ H + <j>>- 8a± (^^1) . (5) 
K dt 2 dt dt k dt \ a dt J 

Eq. ((5]) can be solved with respect to £i{<f>(t)), as 



Combining @ and ©, the scalar potential V(<p(t)) is: 

V(<j>(t)) = ±H(t) a W - 3a(t)H(t)W(t) . (7) 

We now identify t with f((f>) and H with </(£), where / and g are some functions. Such identification has a close 
analogy with the method suggested for the reconstruction of the scalar-tensor theory in Ref . [l9( . Let us consider the 
model where V(<j>) and £i(</>) can be expressed in terms of two functions / and g as 

V(<t>) = ^ff'(/(0)) 2 -^7^- 3 5'(/W))e 9(/(0)) ^) 



1 /■* /'(0 x )e»(/(*O) . 



5' (/W>i)) 



»1 J 



W = f* thtunMe-M+M (lg» (ffa)) + -J- 



(8) 



By choosing V(</>) and £i(0) as (JU), we find the following solution for Eqs.© and (J3]) (compare with (l9|): 

<t> = f-\t) (* = /($) , a = a e^) (H = g'(t)) . (9) 

One can straightforwardly check that the solution ([9]) satisfies the field equation ((4J [l3| ■ 

One should now suggest some realistic ansatze for scale factor and scalar. Our choice below is motivated by the 

possibility to realize the matter dominated stage as well as cosmic acceleration for such scale factor. The choice of 

scalar is dictated by the consistency with scale factor and simplicity condition. Nevertheless, still it remains some 

arbytraryness in our choice so several examples will be discussed. 
Consider now as an example the metric 

c 9{t) = (£T e90t ' /-'(*) = <A)ln^ , (10) 
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where to, go, g\, and <f>o are constants. One further chooses 



,2 2 ffi 

Then U is a constant U — Uo, and 

2 



(11) 



= JL / + 5l e - 0/0o A _ 01-2*/*, _ 3t7() ( + 5l e -0/0o\ jnt/tojoto**'* 



rr /-toe' ■ u / _ \ — / -t 



6w-fy <H w+ t> UJ e ■ (12) 

Eq. Uni) leads to 

H = g + 9 -±. (13) 

Hence, when t is small, the second term in (fT3| dominates and the scale factor behaves as o ~ . Therefore, if 
gi = 2/3, a matter-dominated period, where a scalar may be identified with matter, could be realized. On the 
other hand, when t is large, the first term in (( 13|) dominates and the Hubble rate H becomes constant. Therefore, the 
universe is asymptotically de Sitter space, which is an accelerating universe. The three-year WMAP data are analyzed 
in Ref . [13] , which shows that the combined analysis of WMAP with the supernova Legacy survey (SNLS) constrains 
the dark energy equation of state wde pushing it clearly towards the cosmological constant value. The marginalized 
best fit values of the equation of state parameter at 68% confidence level are given by —1.14 < wde < —0.93. In 
case one takes as a prior that the universe is flat, the combined data gives —1.06 < w^e < —0.90. As in our model 
the universe goes asymptotically to de Sitter space, we find wde — * — 1- Therefore, it can easily accommodate these 
values of wde- For example, if go ~ 40, wde — —0.98. 

In the limit Uo — ► 0, the Gauss-Bonnet term in |T]) vanishes and the action ([I]) reduces into that of the usual scalar 
tensor theory with potential 

V{4>) = \ (go + 2±e-+'*>) 2 - -%e~^° , (14) 

K \ to / K Eq 

which reproduces the result in (2lj . 

Let us now consider a second example. In Einstein gravity with a cosmological constant and matter characterized 
by the EOS parameter w, the FRW equation has the following form: 

Atf 2 = P0 a-^ + . (15) 
Here I is the length parameter coming from the cosmological constant. The solution of (| 1 5(1 is given by 
a = a e 9(t) , 

- d^'-M^ <'-'•»)) ' (16) 

Here to is a constant of the integration and 



o 



2 _ 1 2,2 „ -3(1 



-k I p a . (17) 

We now reconstruct the scalar-Gauss-Bonnet gravity model reproducing (| 1 6|) . If a function g{t) is given by (|16p and 
/(</>) is given by 




/W=*o-^ 7T ^lntanh(-^^— ,1 ^ Ms, 



3(1 + w) \ 4 

U(<p) in ((HI) becomes a constant again, U = Uo- Then, V(4>) and £(<p) are found to be 



V (4>) = -^cosh^ '-0 - \, 2 „ 2 ; smh 2 -V-A- ! -cj> 



2/(3(l+«))) 




-2/(3(l+u)))-l 



(</») = -—^== /,/,>,<■, ,sir- j -^—^ (1J ) { -_sinh( ^ -o, 1 } (10) 
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which again reproduces the result in [21| in the limit of Uq — > 0. Thus, the scalar-Gauss-Bonnet gravity with the 
scalar potentials under considerations reproduces the exact ACDM cosmology. 
Eq. (JTHJ) shows that when f ~ i , the scale factor behaves as a ~ 

(i-io) 2/(3(1+U,)) - Therefore, if to = 0, the 
matter-dominated period can be reproduced. On the other hand, when t —> oo, a behaves as a ~ e*/', which tells 
us that the universe goes asymptotically to de Sitter space, with wde — » — 1. Therefore, it could be consistent with 
WMAP and also with the combined data. 

The Gauss-Bonnet term is usually induced in low-energy string theory. Hence, the string theories could determine 
the form of the functions V(cj>) and £i (cj>). In fact, various types of string compactification give the exponential 
potentials V(tfi) and the exponential function for £i(0). Then the functions V(<f>) and £i(</>) given in this section, 
if they are rather simple, could be given by such compactified string theories. At present, we do not know the 
theory which could give rather complicated functions V(<j>) and £i(</>) as in (|19p. We cannot, however, exclude such 
models, which might appear in future. Such functions might be related with non-perturbative string effects. In our 
formulation, the universe expansion history dictates the possible forms of V(<f>) and £i (</>), which may be obtained 
from (some) string compactifications and the non-perturbative effects. 

III. F(G) GRAVITY RECONSTRUCTION FROM THE UNIVERSE EXPANSION HISTORY AND ITS 

STABILITY 

We can extend the formulation in the previous section to F(G) gravity Q, whose action is given by 
S = J d^x^J—g 

The above action can be rewritten by introducing the auxilliary scalar field as 11 
S = I d x \/—g 



(20) 



(21) 



By variation over 0, one obtains 

Q = V'(<f>)+^)G , (22) 
which could be solved with respect to as 

4> = 4>{G) . (23) 

By substituting the expression (f2"3")l into the action (I2ip . we obtain the action of F(G) gravity, with 

F(G) = -V((f>(G))+^^(G))G . (24) 

Note that the action (|2Tj) can also be obtained by dropping the kinetic term of 4> from the action ((T|) . 

Assuming a spatially-flat FRW universe and the scalar field </> to depend only on t, we obtain the field equations 
corresponding to and ©: 

= - 4^ + V{4>) + 24 ff3^iM) , (25) 



dt 



o = i (m + v{<t>) - w^W)) lQHH d Mm l6H Mm (26) 

k V / dt 1 dt dt 
Combining the above equations, one gets 

K dt 2 dt dt k dt V a dt J 



(27) 



which can be solved with respect to £i (</>(£)) as 

Combining (|23|) and l|28p. the expression for V(<j)(t)) follows: 

VU(t)) = ^H{tf - 3a(t)H(t)W(t) . (29) 
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As there is a freedom of redefinition of the scalar field (j> (which corresponds to the choice f{<j>) = <f> in the last section), 
we may identify t with (p. Hence, we consider the model where V(<f>) and £i(</>) can be expressed in terms of a single 
function g as 

V{4>) = \g'{4>? -3</f»eS«C/(0) , 

K 

1 r<t> e g(4>i) 

U{4>) = ^ [ dfae-^g''^) . (30) 

By choosing V(<f>) and £i(</>) as (|30|) . one can easily find the following solution for Eqs. ([25|) and (|26|) : 

a = a e 9 W (H = g'(t)) . (31) 

After that, one can reconstruct F(G) gravity in a way very similar to the scalar-Gauss-Bonnet theory in the last 
section. 

Although the above formulation is very similar to that in the scalar-Gauss-Bonnet theory, there could be some 
difference: in scalar-GB gravity, since the scalar field has a kinetic term, it can propagate and there can be the 
possibility to generate an extra force besides the Newtonian one. On the other hand, F(G) gravity has no kinetic 
term for the scalar field and an extra force cannot be generated. In fact, one can consider the perturbation around 
the de Sitter background (compare with section five), by writing the metric as g^ u = g^)^ + h^. Here, the Riemann 
tensor in the de Sitter background is given by 

R{0)nvpa = Hi (s(o)ppS(o)i/<7 _ 9(o)iJ.a9(o)up) ■ (32) 

The flat background corresponds to the limit Hq — ► 0. For simplicity, if we choose the gauge conditions g^h^ u = 
V/^-j/fyi, = 0, we find from the equation of motion without the energy-momentum tensor, 

= (V 2 V, - 2Hl hp,) . (33) 

Since the contribution form the Gauss-Bonnet term does not appear except in the length parameter I/Hq of the de 
Sitter space, the only propagating mode should be the graviton in the F(G) gravity. 

Let us now investigate the stability of cosmological solutions in F(G) gravity. The stability study of such theory 
on a de Sitter background will be presented in seventh section where constraints to F(G) from stability condition will 
be defined. For V{<j>) and £i(<£) given by flU, Eqs.(22]) and (25|) yield 

H = 2 + ^ ( 9 "(<P) + g<m , *= ^ + ^e-M (ff 2 - g'W) . (34) 
We now consider the perturbation of solution (f3"Tj) , that is, 

H = g'(t) , <t> = t , (35) 
as 

H = g'(t) + SH , 4> = t + 6<j>. (36) 
Here it is asuumed that 8 H and 6(f) only depend on time t. From (f34|) . we obtain 

Jt ( H ) = M ( % 

M _, 9"'{t) + 4g W'(t) + 



g'(t) ^ k> \ YW 



9"{t) 



If the two real parts of both eigenvalues arc cither negative or vanish, the system turns out to be stable under the 
perturbation. If A is the eigenvalue, the corresponding eigenvalue equations are written as 
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whose determinant is given by 



, g"(t) 2g"(t)e-<>U \ 2 /3 ^A\„ 



The solutions of (l38l) are 



^.'(-^-VM-*^^). (40, 



2 V ff'W 

When Z? < 0, the eigenvalues are complex and, in order that the real parts are negative (as required for stability), we 
have 

A(t)s _m_ V(i) _w)|^ <0 . (41) 

g'(t) k 
When D > 0, both eigenvalues are negative and the system is stable if and only if 

9"(t) „s 2g"(t)e-9W /s / 3 2e~»W\ ,„, , , . 



One may investigate the stability of the cosmological solution ([13)1 : 

= J dtH = g t + g x In 1 . (43) 

Here to is a constant of the integration. Let us consider the case that g\ = 2/3, which corresponds to dust. In the 
early universe, where t — s- 0, we find 

A(t) ~ ^°-t- 8 / 3 > . (44) 
3fv 

Hence, the solution is unstable, which is the typical property of the early-time universe. On the other hand, in the 
late universe, where t —* oo, one finds 

A(t) ~ -4g < , B(t) ~ -4a > . (45) 

Thus, the solution is stable and the future universe does not enter into any Big Rip-like era(22j. 

We also investigate the stability of the ACDM cosmology corresponding to (|16ll . It is convenient to consider the 
case that w — 0, which corresponds to dust. Hence, in the early universe, where t — > to, we find 



2/3 

Therefore, the solution is unstable. On the other hand, in the late universe, it follows that 



Aff ) (i-tor 8/3 >0. (46) 



A ~ ~ < , B(t) ~ 51 e - 3( *-* o)/i > . (47) 
I V 

Thus, the solution is stable again. 

In general, one can consider the case when g(t) behaves asymptotically (that is, in either the early or the late 
universe) as 

g(t) = h In i . (48) 
to 

The case ho — 2/3 corresponds to dust dominated universe. For simplicity, we only consider the case that ho > 
(non-phantom case). Then, we find 

i r 1 - 4/i 2h & h , . . / 3 2 ( t \ ~ ho \ 2h , , 



8 



When < ho < 1/4, it always follows that A(t) > 0; therefore, the solution is unstable. In the early universe, where 
t — ► 0, A(t) behaves as 

As a consequence the solution is always unstable in the early universe. The early universe instability indicates that 
the simple solution (|4"3"]l is unstable and the inflationary epoch should end as is widely expected. 
On the other hand, in the late universe, where t — ► oo, we find 

A^-i-i^, B(t)~!>0. (51) 

Then as long as ho > 1/4, the solution should be stable. 

We have thus presented several examples of the reconstruction program of modified Gauss-Bonnet gravity from the 
universe expansion history which include the epochs of early-time inflation, the matter dominated era, the deceleration- 
acceleration transition, and the acceleration epoch. The stability conditions of such ACDM-like cosmology have been 
investigated, showing that the late universe can indeed be stable, without exhibiting any type of Big Rip behavior 
(for their classification, see [23|). In a similar way, one can consider dark energy cosmology in other versions of the 



hi 



modified Gauss-Bonnet theory [7], [12|, [24( and relate them with the radiation and matter-dominated epochs. 



IV. RECONSTRUCTION OF MODIFIED GAUSS-BONNET GRAVITY WITH MATTER 

It is not difficult to extend the above formulation in scalar-Gauss-Bonnet gravity and F(G) gravity to include 
several matter terms with constant equation of state (EoS) parameters Wi = Pi/pi- Here pi and Pi are the energy 
density and pressure of the i-th matter term. We now include several kinds of matter, like radiation, baryons, cold 
dark matter, etc. Then, instead of ^ and ^ or P5|) and (p?()j) . the corresponding FRW equations are given as 

= -^ 2 + ^ 2 + ^) + E^ + 24 ^^ M ' ( 52 ) 

i 

«- 1 (26 + shA + 1? - v{4>) + y Pi - sfl^ftffi*)) - imii d ^m - leg^yW) . (53) 



n 2 V IT vvv dt 2 dt dt 

Here 77 — 1 corresponds to scalar-Gauss-Bonnet gravity and 77 = to F(G) gravity in the form (|2 1 f) . The energy 
conservation law 

Pi + 3H{p t + Pi )=0, (54) 
gives 

Pi - Pl oa-* {1+Wi) , (55) 
with a constant piQ. Instead of l[8"]). one should consider the model with 

3 ,fttA\\ 2 1_ 

WW 



V&) - -3 , (/(0)) 2 -^T2-3 5 '(/(0))e^))C/ m (^), 



1 /'(0i)e 9(/Wl)) f/ m (^i) 



1_ g'(f(M) 2 



C7 m (0) = j ' dfaf'Me-'UW l^g" (/(^)) + + ^(1 + ^aoe^-O^j . (56) 

In this way we re-obtain the solution ([9]), also in the case when matter is included. This expression is different from 
([5]) due to the last terms in U m , that is, the terms proportional to pio- This term appears since we include matter. 
In |(5]). even without matter, the transition from the matter dominated phase to acceleration can occur. Therefore 
in order that the transition could occur, we need not alway the matters themselves. In the real universe, of course, 
there are matters as in the model in this section. If we consider the model where V{4>) and £i(</>) are given in |56]) 
but matters are not included, the matter dominated phase would not appear. Then in the model (|56|) with matter, 
the matter dominated phase to the acceleration phase occurs not only due to pure Gauss-Bonnet term effects but due 
to combined effect of Gauss-Bonnet term and matter presence. 
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In case of F(G) gravity with 77 = 0, we can choose f(<f>) = 4>. Then, expressions (|56|) can be simplified: 
V(cj>) = l 5 '(0) 2 -3 5 '(0)e^)(7 m (0) , 

^ = 8 J #1 9 W ' 



It is instructive to consider an example of F(G) gravity with only dust as matter, which could be baryons and dark 
matter with w — 0, and g given by 



g(cf>) = I In (sinh (Ctf)) , C ee ~ ^ ^f- . (58) 

This (?(</>) corresponds to (|16p with only dust: u> = 0. Here g(4>) is written in a little bit different way as will soon be 
seen. Eq. (|58f shows that U is a constant U — Uo, and 

2C 2 

VU) = — coth 2 (C0) -3C*C/ocoth(C0)sinh 2/3 (C0) , 

6 = ^ / d0isinh" 4/3 (C(/))cosh 2 (C0) . (59) 



Eq. (|58p indicates that, when <j> = i is small, g(</>) behaves as <?(0) ~ (2/3) In and, therefore, the Hubble rate 
behaves as H(t) = g'(t) ~ (2/3) /t, which surely reproduces the matter dominated phase. On the other hand, when 
(j) = t is large, <?(</>) behaves as 5 ~ (2/3)(C0), that is, if ~ 2C/3 and the universe asymptotically goes to de Sitter 
space. Therefore, the model given by (|59[) with matter shows the transition from the matter dominated phase to the 
acceleration universe, which is asymptotically de Sitter space. In fact, by comparing (|58|) with (|16[) . one can identify 



« = 1, 0="-^. (60) 

We now check if the transition from the matter dominated phase to the acceleration phase could occur or not in 
this model without matter. It can be shown that, if one assumes the existence of the matter dominated phase, a 
contradiction appears. Note that in this model without matter, there is no reason to identify <f> with t. Eq. (|22p gives 

Q = _ 4C 2 coBh(C^) B% ginh -4/3 (c , 0) / 3 _ 2 cosh 2 (C<t>)) + ^ smh- 4 / 3 (C<£) cosh 2 (C<£) G . (61) 
k 2 sinh 3 (C0) V ; 8 

In the matter dominated phase, the curvature R and, therefore, the Gauss-Bonnet invariant should be large. Then 
Eq. (|6"Tj) tells us that <j> should be small in the matter dominated phase, and we find 

G ~ ^g-r 5/3 • (62) 

On the other hand, in the matter dominated phase, the Hubble rate H(t) behaves as H ~ (2/3) /t and, therefore, 

G = 24 (H 2 H + H 4 ) . (63) 
Comparing ^ with $2$, it follows that (j) ~ i 12 / 5 . Hence, 

VW-^-t^, a ~-^3^ 3 ~^ 5 . (64) 

This behavior is in conflict wit h (l2"5"|) . The first equation (|23|) may be written as the usual FRW equation with an 
inhomogeneous EoS ideal fluid [25j (for a particular example of such an inhomogeneous EoS ideal fluid interpreted as 
time-dependent bulk viscosity, see [261 ]) 

> 2 = P G , PG ^V W + 2^ d -^pl. (65) 
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Eqs. (|63f and (|64|) show that pq behaves as i~ 24 / 5 , but H 2 behaves as t~ 2 . Therefore, there is a discrepancy between 
the power of t on both sides of (|65[) . Thus, the matter dominated phase can in no way be realized without matter, and 
consequently the model (|59|) does not generate the transition from the matter dominated phase to the acceleration 
phase without matter. 

As the EoS parameter of the dark energy in the present universe is almost — 1 , the universe could approach the de 
Sitter space asymptotically. Let us consider the form of the action when the universe does become de Sitter space. 
We now assume 

g(t) = H t , f{4>) = M ■ (66) 

Here Hq and /o are constant. Since the Hubble parameter H is given by H = g = Hq, the universe is in fact de Sitter 
space. For F(G) gravity, one may put /o = 1. It follows that 



U (J,) - 71 c -Hofa4> | y Cf + w i)P*0 a (2+3w t )H .f a <t> , tj 



2H f 2 ^ (2 + 3w )H a 

V{(j)) = ^$ + 4 _ y 3(l+»«i)ftoao e3 (i +mi )floW _ 3 Ho U e H °M , 

JO ■ 2 + 

4lW - 16 J ff3/ + Z.24(l + Wj )(2 + 3 Wi )^o 4 8^o 3 1 j 

Here [/o and £o are integration constants. Effectively £o = 0, since a constant times the GB term is a total derivative. 
Even in the presence of matter, the de Sitter universe could still be realized by adding an exponential potential and 
the GB coupling. 

It is interesting to consider the effect of matter to the corresponding cosmology P^|) . This could be done by adding 
the part compensating the contribution from matter to U(<ft), which is chosen to be a constant in the model (|16p . as 



u m (4>) = Uo + Y^W) 



fW ( f3(l+ W ) \ 1 -(4+3uii)/3(l+u0 



Ui{4>) = (l+Wi)p i0 ao J dtiasmhl v % ' (t-t ) ) } (68) 
Then V{4>) and flTSJ) are modified as 



VM = 4 cosh* ?VW^jA 3(1 + .) 
-j \Uo + J2 U i(<P)) cosh 



6(0) 




-2/(3(1+™)) 



8^3(1 + w) 

X cosh^ 2 | K V»v- ■ -/ ^ | I _^ sinh I "V-v- ; -v ^ 1 I ,,,,), 



-2/(3(l+iu))-l 

/ . - , '->/ l i . ,- \ V II I / , - ' '-> i 'I . ,. \ V II 



Therefore, there exists a scalar-GB gravity which shows the transition from the matter dominated phase to the 
acceleration phase, even in the presence of matter. This theory could be regarded as an extension of (|58[) with only 
dust (baryon and/or dark matter). Thus, in the model (|68[) without matter, the matter dominated phase cannot be 
realized and therefore the transition from the matter dominated phase to the acceleration era does not occur. 
Let us investigate the asymptotic behavior of Ui(4>) in ([68]) . When t = f(<f>) — » to, one gets 

w,— « ^Mi^ (/w _ (o) l'-'«'---~ . (70) 



a{-l + 3(w-Wi)} \ 21 

When the model containing the matter dominated phase is considered, one can find w = 0. Since usually Wi > 0, the 
power [—1 + 3(w — Wi)] /3(1 + w) could be always negative. Hence, Ui(<j>) could be singular when t — f{<j>) — ► to- This 
is not something extraordinary. For example, the scale factor of the FRW universe with matter, whose EoS parameter 

2 

is a constant w, behaves as a ~ t 3 ( 1 +«") , which is singular at t = 0. 
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On the other hand, when t = f(<fi) — > oo, we find 

/n>\ -(4+3uii)/3(l+u>) 9/ / 4 + 37/;- \ 

w> 4rk^(-^r 1(/W) -'° ) ) ' (71) 

Here Uip is an integration constant. When t = /(</>) — ► oo, the second term vanishes very rapidly. We may choose 
UiF = 0. Thus, in the case Uo ^ 0, the corrections generated by adding matter could be neglected. These corrections 
become important only when t — f(<p) — > to, which may correspond to the early universe. 

Our study shows the existence of a big class of modified Gauss-Bonnet gravity models with matter where the 
cosmological sequence of matter dominance, deceleration- acceleration transition and cosmic acceleration occur only 
in the presence of matter. 

V. A COMPENSATING DARK ENERGY 

Some versions of scalar-Gauss-Bonnet or /(G) gravity (with specific potentials) do not have a matter dominated 
stage as a solution, not even in the presence of matter. Similarly, some of them have a stable matter dominated era 
as a solution, hence, no transition to an acceleration epoch occurs. It is remarkable that even in this case, a realistic 
cosmology may emerge at the price of introducing a compensa ting dark energy. Such scenario has been proposed 
in Ref. [lg], based on the example of modified gravity of Ref. [27|. Let us introduce a compensating dark energy, 
which could be an ideal fluid and may help to realize the matter dominated and deceleration-acceleration transition 
phases in modified GB gravity. One may question the necessity of such ideal fluid, saying, that already the Einstein 
gravity with negative pressure ideal fluid may describe the acceleration of the universe. The below scenario is more 
sophisticated: it is supposed that the acceleration is produced by scalar-GB interaction and compensating dark energy 
is negligible in the late epoch. (Hence, without GB term no cosmic acceleration occurs). The role of compensating 
dark energy is relevant only in matter dominated era: it helps in the realization of matter dominance. 

When the energy density p c and p c of the compensating dark energy are added, the FRW equations become 

= - '^H 2 + l4> 2 + V{<j>) + 24H M f} t)) + Pc + Pm , (72) 



K~ 



2 v ' dt 

j2 



= \(2H + <IhA +V -y(0)-8g 2 ' - 16gg ^#M - Ifig* +Pc+Pm . (73) 

k V / 2 dt dt dt 

Here, p m and p m are contributions from matter, which may also include dark matter. In the matter dominated phase, 
the contribution from the scalar-GB terms could correspondingly be canceled by that coming from the compensating 
dark energy. In other words, in the matter dominated era the following conditions hold 

Pc ~ -{^ 2 + ^) + 24i/3^i|M| , (74 ) 

Pe ~ - to - vw - zh^w)) _ imk <Mm _ x ^mm \ . (78) 

^ 2 dt 2, dt dt J 

If matter is dust, as baryons or cold dark matter, the Gauss-Bonnet invariant G behaves as (|(kJ]) . 

For simplicity, we will first consider F(G) gravity with r\ = 0. When H = h/t, G is given by G — 24h 3 (h — l)/£ 4 . 
Then G changes its sign when h = 1. In the decelerating universe, h < — 1, but in the accelerating one the sign of h 
changes: h > — 1. When the decelerating universe turns to the accelerating phase, G changes its sign. In this case, 
solving (|2"2")l . if we do not choose V(</>) and £i properly, there is no solution. In the case H = h/t, Eqs. (f5D|) yield 

t r/ / \ 3fe 2 (fe~l) ^ 

VW= (ft + l)«V 1 ^ ) = 8h(h+iW (76) 

We now consider a new model by replacing h in (|76[) with h(cj>), adiabatically depending on 0, so that one can neglect 
all derivatives h'(4>), h"(4>), • • •: 



ih^fW) - 1) 

(h(4>) + i)k 2 4> 2 ' ^ v-r/ 8/1(0) + i)k : 

Using (|22|) and neglecting /i'(0), h"{(j>), • • •, we find, as expected, 

^ 2 = 24fe(^) 3 (l-feW) _ (?8) 
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As the adiabatic approximation is used, our identification <f> = t can again be introduced. For the matter dominated 
era, by using (f74"|) and ([75]) , we find 



3h(t) 2 h(3h(t) - 2) 

k 2 4>- 



which gives the EoS parameter 

p c 2 — 3h , „, 

Hence, when /i — > 2/3, which corresponds to dust, w c —> as expected. We should note that the energy density of 
the compensating dark energy p c is now negative, which might be possible if there is a negative cosmological constant 
(for instance, the one produced by anti-deSitter space), which effectively shifts the energy a by negative constant. 
We now consider the scalar-Gauss-Bonnet theory (with n = 1) of Ref. [5[, where 

V = V e~% , &(<W=£oe^. (81) 
Without matter, the theory with those potentials (|81[) can be solved [|[ 
H = if , = 0o In £ , (when ft„ > 0) , 

H = -£ti , = O In ^ , (when ft < 0) . (82) 



Here ho and ii are related with Vq and £o, as 

^.--l fop <!-«•> 

Av (1 + Alo) 

48Co^ 6 / 



K 2 (1 + fto) 



fto-^r • (83) 



Since the above ho is a constant, the EoS parameter w is also a constant. Therefore, in this model without matter 
and compensating dark energy, the transition from the matter dominated phase to the acceleration phase could not 
be generated. 

In the following, it is enough to consider the case with h a > only. Let us now add the contribution from the 
compensating dark energy and usual matter, and assume H(t) = h(t)/t with a slowly varying function h(t). The 
time dependence of 4> can be determined by solving the scalar field equation Neglecting h'((f>), h"(4>), ■ • • in the 
adiabatic regime, we get 

= o l n -^- Y . (84) 
ti(t) 

Here t\(t) is found from the following equation 



= (1 - 3h{t)) <& + 2V h(t) 2 - 48g °y )3 (h(t) - 1) . (85) 



n 

Hence, in the matter dominated phase, Eqs. (|74|) and ([75]) are 

Pc ~ -{f + V + m a h(tf}± , (86) 

Pc ~ - | y -V + moHt) 2 (1 - 2fe(t))| ± , (87) 
which give the effective EoS parameter 

f - V + moHt) 2 (1 - 2h(t)) 
w — — 

f + V +48£oh(t) 3 



(88) 



Thus, scalar-Gauss-Bonnet gravity with potentials (|81[) and without matter and/or compensating dark energy does 
not contain the transition from the matter dominated phase to the acceleration phase. Nevertheless, adding the 
matter and the compensating dark energy (|86p . the matter dominated phase could be realized even for the model 
(|8"Tj) . similarly to other classes of modified gravity [l6j]. As expressed in ([74]) . the compensating dark dark energy 
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cancels the contributions from the scalar-Gauss-Bonnet term in the early universe and, therefore, the system effectively 
reduces to Einstein gravity coupled with matter, where the matter dominated universe can in fact be generated. In 
the late universe, the contributions from matter and the compensating dark energy become small and the accelerated 
expanding universe naturally takes over. A similar scenario can be devised to improve the emergence of the matter 
dominance stage in other versions of modified gravity and scalar-tensor theories (for a recent study of accelerated 
cosmologies in scalar-tensor theories, see [28[ and the references therein). 

VI. DISCUSSION 

In summary, we have tried to show in this paper that string-inspired scalar-Gauss-Bonnet gravity and modified 
Gauss-Bonnet gravity are indeed interesting alternatives for dark energy. These theories are in fact closely related 
and, what is also important, may have a stringy origin. Here, the reconstruction program from the universe expansion 
history for those theories has been carried out successfully. Several explicit examples of the same (with some specific 
potentials) have been presented where the cosmological sequence of the matter dominance, deceleration-acceleration 
transition and cosmic acceleration occurs very naturally. Moreover, the accelerated universe can be asymptotically 
de Sitter or it may correspond to an exact ACDM cosmology. The study of perturbations around the cosmological 
solutions above has been performed too. There is no problem to include usual matter with specific equation of state 
into such consideration. In that case, one can construct versions of the Gauss-Bonnet gravities above where a matter 
dominance period and a deceleration-acceleration transition occur only in the presence of matter. It is also remarkable 
that, even in the case when such intermediate universe is not a solution of some modified Gauss-Bonnet gravity, it 
can actually be made so, at the price of introducing some compensating dark energy. 

Since the Gauss-Bonnet term is usually induced in low-energy string theory, the functions V{<j>) and £i(0) could con- 
tain the information about the string compactification and/or stringy non-perturbative effects. By the reconstruction 
program in this paper, we suggested the possible forms of V(<p) and £1 (</>), which are cosmologically viable. 

Due to the fundamental role of the de Sitter space which appears in our scenario as the final state of the universe, 
special attention has been paid to such space. In particular, the one-loop effective action of F(G) gravity has been 
found on the de Sitter background. This effective action was then used to derive stability criteria for the modified 
Gauss-Bonnet gravity theory. Some numerical examples showing explicit versions of modified Gauss-Bonnet gravity 
with a stable de Sitter vacuum have been presented. 

The successful reconstruction of string-inspired, Gauss-Bonnet gravity from the universe expansion history, per- 
formed in the present work, shows that it actually represents a reasonable gravitational alternative for dark energy. 
Having in mind the promising results obtained in the comparison of such a theory with observational data (see, for 
example, @), it becomes clear that it deserves careful attention. Moreover, the theory successfully passes the check 
of the three-year WMAP observational data. Needless to say, as with any other alternative to General Relativity at 
the current and/or future universe, additional accurate checks regarding the solar system tests should still be carried 
out. Nevertheless, even if some problems could be encountered, there will always remain a reasonable chance that the 
situation improves by taking into account higher-order string corrections, as has been wisely indicated in Ref. [i"3| . 
The introduction of such higher-order string corrections in the above scenario will be discussed elsewhere. 
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APPENDIX A: ONE-LOOP EFFECTIVE ACTION IN F(G) GRAVITY ON DE SITTER SPACE 

The important issue in any fundamental gravitational theory is stability issue, as it indicates if some of the highly 
symmetric spaces (flat, de Sitter or Anti-de Sitter) could be the ground state. Having in mind that de Sitter space 
appears as classical solution of the above theory in the early as well as in the late universe we study stability of GB 
modified gravity in de Sitter space. 

In order to study the stability of F(G) gravity, one possibility is to calculate the one-loop effective action on the 
corresponding background. It will be here shown that, in the calculation of the one-loop effective action in higher- 
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derivative gravity a multiplicative anomaly appears naturally. Thus, before starting the calculation we first review 
this concept in a general situation. The usual definition of the determinant of an operator, A, is done through its zeta 
function [29| Ca(s) = Yliel ^T* = * r A ~ s , where analytic continuation in s to cover the domain of the complex plain 
to the left of the abscissa of convergence s = sq is naturally assumed. Recall that sq — m/n, namely it is given by the 
dimension of the working space, m — dim M, divided by the order of the operator, n — ord A. For an elliptic operator 
of positive order on a compact manifold, satisfying the usual spectral condition (e.g. the presence of an Agmon- 
Nirenberg cut in the spectrum), the zeta function exists and has very nice general properties, as being meromorphic 
with possible single poles on integer and specific fractional positions Sk = (n — k)/m, k = 0, 1, 2, . . . , n — 1, n + 1, . . . 
of the negative real axis only [3(| • The determinant of A is then given by [3l| 



det c A = exp[-Ci(0)]. 



(Al) 



This definition only depends on the homotopy class of the spectral cut. 

Now, given the operators A, B and AB, even if (a, Cb and (ab exist, it turns out that, in general, det^(AB) ^ 
dct^A dct^B. The multiplicative, or noncorrrmutative, or determinant anomaly (also called defect of the determinant) 
is defined as: 



S(A, B) = In 



det c (AB) 



det c A dett B 



C^(Q) + Ca(o) + Cb(o)- 



(A2) 



There is a useful formula due to Wodzicki for the multiplicative anomaly [32, 33]. In terms of Wodzicki's residue: res 
A = 2 Res s= o tr(AA~ s ), A being the Laplacian operator, or equivalently, in local form using the symbol expansion, 
res A = Js« M tr a_„(:r, £) ci£, with S*M C T*M the co-sphere bundle on M, the multiplicative anomaly can be 
obtained as 



8{A,B) 



res{[lncr(A,S)] 2 } 
2 ord A ord _B(ord A + ord B) ' 



a(A,B) :=A ordB B -ordA 



(A3) 



Several implications of this multiplicative anomaly have been extensively discussed in the literature [34|. Here we 
consider a very different situation, never met before, namely the case when one of the two operators slowly vanishes 
owing to the change of a certain parameter, a situation that may be not too unfrcqucnt. The main issue here is to 
check if, in fact, the shrinking operator could leave some imprint on the surviving operator, through the multiplicative 
anomaly term [341 ]. To simplify notation we will drop the C from the determinants hereafter. 

Such anomaly appears basically in higher derivative gravities, like the R 2 theory, when one calculates the one- 
loop effective action in some background. This will precisely be case under consideration here. The anomaly may 
persist even after one of the operators shrinks, say as a consequence of time evolution or by the action of some other 
parameter, as we are now going to see. In fact, as an example, let us consider the uniparametric family of operators 
B e , where e is a real parameter that will eventually go to zero, e.g., adiabatically. For further simplicity, let us just 
take that B e — eB, an overall constant factor. It is easy to see from its definition that the anomaly S(A, B e ) does not 
actually depend on e, in fact 



S(A,B e )=6(A,B). 



(A4) 



And, again from the definition of the anomaly (|A2[) , it turns out that in the limit when the B e operator adiabatically 
disappears, a contribution may remain, in the way that 



det(AB e ) 



dct B f 



detA-e 5(A ^ 



(A5) 



This is again an implication of the presence of the multiplicative anomaly and, on its turn, of the definition of the 
zeta determinant. It contributes, generically, an additional term to the determinant of the resulting operator after 
taking the limit. 

Now to the physics. Corresponding to the scalar version of F(G) gravity, we have the Euclidean action 



S E [g] 



1 



R 



eg ij di4>dj4> 



F(4>)+F'(4>)(G-4>) 



(A6) 



which is equivalent to the modified Gauss-Bonnet gravity as investigated in [14| . In (|A6|1 (j> is a scalar field which "on 
shell" becomes equal to the Gauss-Bonnet invariant G, gij is the metric, R the related scalar curvature and finally 
F ((/)) is an arbitrary smooth potential. An interesting observation is that the multiplicative anomaly commutes with 
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the limit of the small parameter e — > 0. We will confirm that this is the case below. It may persist after the limit is 
enforced. 

As it as been shown in [3] , to ensure the existence of de Sitter solutions the function F has to satisfy the condition 

GF'(G)-F(G) = j, (A7) 

where the quantities G and R are respectively the Gauss-Bonnet invariant and the scalar curvature related to the de 
Sitter constant solution. 

Now we are going to consider small fluctuations around the constant curvature solution gij, associated with de 
Sitter solution, so we put 

fjij = 9i 3 + hij , h = g v hij , 4> = G + (f> , (A8) 



and we develop the action (|A6[) around the constant curvature solution up to second order in and <f>. As usual all 
tensor indices are lowered and risen by means of . After a straightforward calculation, taking into account that we 
are dealing with a maximally symmetric space, for the quadratic part and disregarding total derivatives we obtain 

- -fohfj R h % , IiR 2 h 2 tj hij A hij faRhijAhij 

C2 = — 6~ + - 6~~ + + 6 

, foRtf , R 2 tf 31 /i R 3 & 2 , / &A& 



8 16 288 2 

R&Ati 17/xi? 2 &A& fx R^AA^i 



4 36 6 

foh 2 f 2 <t> 2 fih 2 R 2 hh^pR 2 ShAh 

16 2 48 6 32 

hhhRAh e<j>A(j> j 2 hRA<j) 

32^ + 2 A 
hRAa f 1 hR 2 Aa f 2 <f>R 2 Aa 

16 16 + 12 

foRaAa R 2 a Act 17/ii? 3 crAcr 

16 32 + 288 

3/iAAcr 3 /i h R A A a f 2 <t>RAAa 

16 16 + I 

3/ocrAAcr RaAAa 3f 1 R 2 aAAa 

16 8^ + 16 

3cr A A Act fiRa AAAa 

32 + 32 ' 



(A9) 



where ff. — F^ k '(G) (k — 0, 1, 2) is the derivative of the function F(<j>) evaluated at G, while hij, a and h are the 
irreducible components of the symmetric tensor field hij . They are related by 

h %3 = h^ + V^j + + ViVjcr + - g l3 (h-Aa), (A10) 
and satisfy the conditions 

V fe £ fc = 0, Vih ij =0, g ij hu = 0. (All) 

Here and A = £f* J 'V<Vj are respectively the covariant derivative and the Laplace-Beltrami operators related to 
the de Sitter metric gij . It has to be remarked that Eq. (|A9[) coincide with the analog expression written in [l4| when 
one puts (j> = G and e = (see Sec. IV in the cited paper). 

As it is well known, invariance under diffeormorphisms renders the operator in the (h, a) sector not invertible. One 
needs a gauge fixing term and a corresponding ghost compensating term. We consider the class of gauge conditions 

Xk = Vjhjk - v fc h > 

parametrized by the real parameter p. As gauge fixing we choose the standard expression 

£ g f = lx i G ij x i , Gij=a gi j (A12) 
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and so the corresponding ghost Lagrangian reads [3f 



5 y 3 , 



(A13) 



where Ck and Bk are the ghost and anti-ghost vector fields, respectively, while S x h is the variation of the gauge 
condition due to an infinitesimal gauge transformation of the field. It reads 



5 hi 



= gij A + Ri 



Neglecting total derivatives one has 



R 



£ gh =aB« [ A + — ) C k 



(A14) 



(A15) 



and finally in irreducible components we obtain 

2 

r (a 



c 



R 



R 



4/ ^ + iH A + i Ufl 



9 



16 



p " hAh-—a [A 
16 



i? 



A a 



(A16) 



Ac 



where ghost irreducible components are defined by 



Cfe — Cfe + VfeC, 
S fe = B k + V k b, 



v k c k = o. 



V k B k = . 



(A17) 



(A18) 



In order to compute the one-loop contributions to the effective action one has to consider the path integral for the 
bilinear part 



C — C2+ Cgf + Cgh 



(A19) 



of the total Lagrangian and take into account the Jacobian due to the change of variables with respect to the original 
ones. In this way, one gets [H, [36| 



Z<U = (det G l3 y 1/2 J D[h^]D[C k ]D[B k ] exp (- J d A x^C 



1/2 



(detG y ) 1/2 det Jf 1 det J, 



£»[/i]£>[%]£)[e]£>H J D[C fe ] J D[B fe ]£»[c] J D[6] exp - / d^x^C) , 



(A20) 



where the determinant of the operator Gij for our choice is trivial, while J\ and J2 are the Jacobians due to the 
change of variables in the ghost and tensor sectors respectively. They read (35j 



Jx = A, 



A, fj f-A -~] A„. 



(A21) 



A and A 1 being the Laplacians acting on scalar and vector fields. 

Now, a straightforward computation (disregarding the multiplicative anomaly ([13])) leads to the following off-shell 
one-loop contribution to the "partition function" 



- r(1) = = 



det 



" Al -4 



x det 



-A 2 - 



1/2 , ( a R 
det f -A -- 

R{9f a +4R- 12X) 
3(4/o + 3R- AX) 



-1/2 
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det 



9e 21e/ f 2 2 R 2 21eX 



A 3 
^ 



256 256 R 512 256 R 
VTefp 9/ 2 21 / / 2 7efl 7/ 2 2 ^ 3 31eA 21 / 2 X 
256 256 256 i? 256 1536 512 256 R 



A 



17/o/ 2 _ 9e/ i? 7/ 2 fl _ 3ei? 2 5 / 2 2 R 4 _ 31f 2 X 3eflX ' 
256 1024 256 1024 1536 512 512 ! 

9/0/2^ 3/ 2 fl 2 / 2 2 i? 5 3f 2 RX " 1 ' 2 

1024 1024 1536 512 



(A22) 



For convenience we have written this in the Landau gauge corresponding to p = 1 and a = 00 and moreover we have 
put X = fo + R/2 — Gf\. Thus, X = is the "on-shell" condition corresponding to the de Sitter solution. 

As a check one can easily see that Eq. (1A22[) has the correct limit in the Einstein plus cosmological constant case, 
that is when s = 0, fi = 0, /o = —1/2. In fact in such a case the complicated expression in the scalar sector decouples 
and one obtains the well known result reported in Ref. [36J . Furthermore, when e = 0, it can be shown that also 
at one-loop level such a model is equivalent to the F(G) theory developed in [l4l |. which confirms the commutativity 
with the multiplicative anomaly. 

APPENDIX B: STABILITY CRITERIA FOR MODIFIED GAUSS-BONNET GRAVITY 

In this Appendix, applying the results of the previous Appendix we study the stability issue for several specific 
models. In principle, one can study the one-loop effective action for the de Sitter space explicitly in terms of special 
functions [36j |. This happens in some cases, where one is able to express the determinant, say det (— Aq + aiA 2 , — 
a 2 Ao + 03), in terms of a product of more elementary determinants of lower-dimensional operators. It is there that 
the multiplicative anomaly explicitly appears. Fortunately, for the case under discussion we are able to carry out a 
detailed analysis of the stability conditions of the model without the need to go through all this cumbersome process. 
This is what we do here. 

In reference [3a |. the one- loop effective action has been also used in deriving a stability criterion for the class of 
modified models described by the function f(R). Such a criterion of stability has been obtained by imposing the 
absence of negative eingenvalues related to the Laplace like operators appearing in the regularized one-loop effective 
action and it has been independently confirmed in [39l | within a pure classical approach, namely, involving the study 
of cosmological perturbations. 

In the case we are discussing here the situation is more complicated from the technical view point, since the 
functional determinant one has to investigate involves an algebraic polynomial of third order in the Laplace operator. 
But, in principle, a stability criterion for such a class of modified Gauss-Bonnet models could be investigated along 
the same lines, namely by requiring the vanishing of the imaginary part in the one-loop effective action or the absence 
of negative eingenvalues. Note that the equivalent analysis of cosmological perturbations in F(G) gravity is extremely 
complicated [9J. However, such study may be very important cosmologically, since it may prove (or disprove) the usual 
belief that the current (almost) de Sitter dark energy era may be eternal. 

In order to discuss this issue we take the limit e — > and consider the on-shell condition X = 0. In principle the 
operator of the scalar sector in (|A22[) may be written in the factorized form 

A(X) = X 3 + a x X 2 + a 2 X + a 3 = {X - X X ){X - X 2 )(X - X 3 ) , (Bl) 

where now X = — Aq is a non-negative differential operator, Xi are the roots of the third order algebraic equation 
A(X) = and the constant coefficients ai are given by 

f7 18 , 42/o 
ai - R {3 + 7^ + frR* 
R 2 (5 , 14 34/o 

a 3 = -R 3 



3 f 2 R 3 f 2 R 4 
1 , 3 , 9/o 



^3 2/ 2 i? 3 f 2 R 4 , 
The nature of the roots depend on the discriminant, D, which reads 

D _f 3a 2 -ai V | / 9a 1 a2-27a 3 -2 a 3 y (m) 
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Depending on the sign of the discriminant one has the three cases: 

(i) D > 0, then one root is real while the other two form a couple of conjugated complex numbers, say X\ = A3. 

(ii) D = 0, then all the roots are real and at least two of them are equal, say X 2 = X3. 

(hi) D < 0, then the three roots are real and distinct. 

As a consequence, in the first two cases the factorization reads A(X) = (X — Xi)B(X), B(X) being a non negative 
differential operator and then one has stability of the Gauss-Bonnet modified models provided X\ < or, which is 
equivalent, 03 > 0, since X\X 2 Xj, = —03. Then requiring 03 > and D > we arrive at the following sufficient 
conditions for the stability of the model we are considering on the de Sitter background: 

<i.h-R 2 h-3R = u, 
J2 h 

76204800 /o 4 + 136950912 / 3 R + 92 1 61152 f 2 R 2 

+27527040 f R 3 + 3079296 R 4 + 690 1 200 f 3 f 2 R 4 
+9430128 f 2 h R 5 + 4292784 f f 2 R e 

+651024 h R 7 + 171975 f Q 2 f 2 2 R s + 154794 /„ f 2 2 R 9 

+34815 f 2 R 1Q + 1600 f f 2 3 R 12 + 704 / 2 3 R 13 < . (B4) 

The first one is the on-shell condition, while the second and the third ones derive from 03 > and D > respectively. 

When the discriminant is negative all the roots are distinct and the condition 0,3 > does not ensure the operator 
A(X) to be non-negative. Of course, a sufficient condition can be obtained by requiring all the roots to be negative, 
but in such a case, for technical reasons we did not find a reasonable stable model. 

As an example let us consider the choice F(G) — a\fc + /3G 7 , a and 7 being arbitrary dimensionless parameters, 
while f3 has dimensions [mass] 1-27 . This is the only dimensional parameter in the model and for this reason it remains 
free also after imposing the above restrictions. The on shell condition gives 

a*-. £3^33 , ( B5) 

where a, (3, 7 are assumed to have values such that R > 0, since we want a de Sitter solution. The other two conditions 
above fix the possible ranges of a and 7. In particular it can be seen that there are stable solutions for very small, 
positive values of both a and 7, only. 

A careful analysis of the stability conditions (|B4[) leads to the following conclusions. First, the whole analysis can 
be done in terms of /o and f 2 , by the substitution of the first equation, namely 

'■=!(tM- (B6) 

Then, an analysis of the roots of the last inequality, seen as a polynomial of fo at the bound of the second inequality, 
shows that there is a stability region for values of /q around fo — —R/2, that is f 2 ~ —6/R, and f 2 < 9/4i? 3 . This 
is depicted in the figures. All of them are plots of the last inequality in (|B4[) . where the relevant domain for fo is 
identified (as we see there, an additional one exists but it requires very precise fine tuning). Fig. la is obtained for f 2 
on the first border of the first of the inequalities in (IB4|) . which is two-sided, namely 



f 2 < (-9/2i? 3 )(l + 3/o/fl), fo < -R/3, (B7) 
h > (-9/2i? 3 )(l + 3f /R), fo > -R/3. (B8) 

Fig. lb shows the same plot but for f 2 well inside the region delimited by the first of the inequalities: f 2 — 1/R 3 . 
Fig. 2a shows the stability region corresponding to the very particular case of a vanishing second derivative of F(G), 
namely f 2 — 0, while Fig. 2b corresponds to the value f 2 = —1/R 3 , which marks quite closely the end of the first 
stability domain (|B7[) for f 2 . We thus see that for very reasonable values of the first derivatives of the function F(G), 
we are inside the domain of convergence. To identify specific functions having these derivatives, however, is not an 
easy task, as the preceding example shows. 

Similarly, one can study the stability conditions of the de Sitter universe which corresponds to the final era in the 
first of the reconstruction scenarios considered above, for other versions of string-inspired gravity. It is our impression, 
that this study, which recently became quite important due to the possible emergence of a (metastable) de Sitter 
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FIG. 1: Fig. la, left, plot of the lhs of the last inequality in (|B4|) . which we call ^3, where fa lies on the left border of the first 
stability domain (|B7[) . Fig. lb, right, same plot where /2 = 1/R 3 lies inside the region delimited by the first of the inequalities 
in (|B4[I . not far from the right end of the first stability domain 
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FIG. 2: Fig. 2a, left, plot of Is, where /2 corresponds to the case /2 = 0. Fig. 2b, right, same plot where /2 = —1/R 3 lies close 
to the left end of the second stability domain (|B8f) . 



vacuum in string theory, could prove to be less involved than the corresponding analysis of cosmological perturbations 
@ . Indeed, even the classical stability study of higher-derivative gravity on the de Sitter background is already not 
that easy [40j ]. 
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